Damping of filament thread oscillations: effect of the slow continuum 
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ABSTRACT 

Transverse oscillations of small amplitude are commonly seen in high- 
resolution observations of filament threads, i.e. the fine-structures of solar fila- 
ments/prominences, and are typically damped in a few periods. Kink wave modes 
supported by the thread body offer a consistent explanation of these observed oscilla- 
tions. Among the proposed mechanisms to explain the kink mode damping, resonant 
absorption in the Alfven continuum seems to be the most efficient as it produces damp- 
ing times of about 3 periods. However, for a nonzero-/? plasma and typical prominence 
conditions, the kink mode is also resonantly coupled to slow (or cusp) continuum 
modes, which could further reduce the damping time. In this Letter, we explore for 
the first time both analytically and numerically the effect of the slow continuum on the 
damping of transverse thread oscillations. The thread model is composed of a homo- 
geneous and straight cylindrical plasma, an inhomogeneous transitional layer, and the 
homogeneous coronal plasma. We find that the damping of the kink mode due to the 
slow resonance is much less efficient than that due to the Alfven resonance. 

Subject headings: Sun: oscillations — Sun: magnetic fields — Sun: corona — Sun: 
prominences 



1. Introduction 

The existence of the fine- structure of solar filaments/prominences is clearly shown in high- 
resolution observations by both ground-based and satellite-on-board telescopes. These fine-structures, 
here called threads, have a typical width, W, and length, L, in the ranges 0.2 arcsec < W < 
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0.6arcsec and 5 arcsec < L < 20arcsec (|Lin et al.ll2005|) . Filament and prominence threads are 
believed to be thin magnetic flux tubes partially filled with prominence-like plasma and whose 
footpoints are anchored in the photosphere. Obse rvers usually report the presence of transverse os - 
cillations and waves in these fine-structures (e.g. lLin et al. Il2003[l2005[l2007l:lokamoto et alJ2007h. 



which have been interpreted in te rms of magnetohyd r odynamic (MHD) waves (e.g., 



Terradas et al. 



2008). The reader is referred to bliver & Ballesterl d2002h . ballesterl d2006h . and IBanerjee et al. 



(I2007|) for reviews of these events. In addition, a typical property of small-amplitude prominence 
oscillations observed in Doppler series is that they are quickly damped, w ith damping times of the 
order of several periods (|Molowny-Horas et al.ll 19991 : iTerradas et alj|2002l) . Since small- amplitude 
oscillations are of local nature and related to periodic motions of threads, there must exist some 
mechanism which is able to damp thread oscillations in an efficient way. 

The damping of prominence oscillations h as recently been in vestigated in a number of papers 
(this topic has been reviewed by IOliverll2008|) . ISoler et al.1 (I2008|) modeled a filament thread as a 
homogeneous cylinder embedded in an unbounded corona and studied the effect of nonadiabatic 
mechanisms (thermal conduction, r adiative losses, a nd heating) on the wave damping. As in pre- 
vious works (jCarbonell et al.ll2004l : ISoler et al.ll2007l) . they concluded that thermal effects are only 
efficient in damping the slow mode, the fast kink mode being almost undamped. Subsequently, 
Arregui et al.l (I2008|) used a more realistic thread model and included an inhomogeneous transi- 
tion region between the cylindrical thread and the external medium. He neglected plasma pressure 
and adopted the so-called /? = approximation, where (5 is the ratio of the plasma pressure to 
the magnetic pressure. These authors found that the kink mode is resonantly coupled to Alfven 
continuum modes due to the presence of the transverse inhomogeneous layer, and obtained a ratio 
of the damping time to the period td IP ~ 3 for typical prominence conditions. In th e case of 
coronal loops (see, for example, recent reviews by Goossens et al. 2006 : Goossensll2008 ). resonant 
absorption has been extensively investigated as a damping mechanism for the kink mode. 

Although the plasma J3 in solar prominences is pr obably small, i t is de finitely nonzero. Hence, 
if gas pressure is taken into account in the model of I Arregui et al.1 (|2008|) . i.e. the /? case, it 
turns out that the kink mode phase speed, namely ct, is also within the slow (or cusp) continuum 
that extends between the internal, c si , and external, c se , sound speeds, i.e. c si < q. < c se . In this 
case, the kink mode is not only resonantly coupled to Alfven continuum modes but also to slow 
continuum modes. The frequency of the kink mode is both within the Alfven and slow continua 
because of the high density and low temperature of the prominence material in comparison with 
the coronal values. Hence, the slow continuum damping arises as an additional mechanism to 
damp the kink mode in filament threads, and its efficiency needs to be assessed. Although the slow 
resonance has been previousl y investigated in t he context of absor ption of driven MH D waves 
in the solar atmosphere (e.g.. ICadez et al.l 1 19971 : Erdelyi et alJboOll) and sunspots (e.g.. Keppens 
1996 ). the present work studies for the first time the effect of the slow resonance on the damping of 
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the kink mode in filament threads. For coronal loops the kink speed is outside the slow continuum. 
A coronal loop is presumably hotter and denser than its surrounding corona, so the ordering of 
sound, Alfven, and kink speeds is c se < c s; - < v Ai < < v Ae . Therefore, there is no slow resonance 
for the kink mode in coronal loops and the present study does not apply to coronal loops. 

This Letter is organized as follows: § |2]contains a description of the model configuration and 
the mathematical method. The results are presented and discussed in § [3j Finally, our conclusions 
are given in § |U 



2. Model and method 



The model for the equilibrium configuration considered here is equivalent to that of lArregui et al 



(120081) . We model a filament thread as a straight cylinder with prominence-like conditions embed- 
ded in an unbounded corona, with a transverse inhomogeneous layer between both media (see 
Fig.Q]). We use cylindrical coordinates, namely r, <p, and z for the radial, azim uthal, and longitudi- 
nal coordinates. We adopt the density profile by iRuderman & Roberts! (|2002l) which only depends 
on the radial direction, 



with 



Po O) 



P,(r) = - 




r<a- 1/2, 
1/2 < r <a + l/2, 
r>a + 1/2, 



sin 



71 ( 



a) 



(1) 



(2) 



Pi I \ Pi 

In these expressions, p, is the internal density, p e is the external (coronal) density, a is the tube 
mean radius, and / is the transitional layer width. The limits l/a = and l/a = 2 correspond to a 
homogeneous tube and a fully inhomogeneous tube, respectively. We use the following densities: 
Pi = 5 x 10" 11 kg m" 3 and p e = 2.5 x 10" 13 kg m" 3 . Therefore, the density contrast between 
the internal and external plasma is Pi/p e = 200. The plasma temperature is related to the density 
through the usual ideal gas equation. We consider T, = 8000 K and T e = 10 6 K for the internal 
and external temperatures, respectively. The magnetic field is taken homogeneous and orientated 
along the z-direction, B = B z, with B Q = 5 G everywhere. With these conditions,/? ~ 0.04. 



2.1. Analytical dispersion relation 



We consider the linear, ideal MHD equations for the B + case. Since (p and z are ignor- 
able coordinates, the perturbed quantities are put proportional to exp (icot + imp - ik z z), where co is 
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Fig. 1. — Sketch of the model configuration. 



the frequency, and m and k z are the azimuthal and longitudinal wavenumbers, respectively. In the 
absence of a transitional layer, i.e. l/a = 0, the dispersion relation of magnetoacoustic waves is ob- 
tained by imposing the continuity of t he radial displacement, £ r , and the total pressure perturbation, 
Pi, at r = a. This dispersion relation (Edwin & Robertsll 19831) is T) m (co, k z ) = 0, with 



D m (co, k,) = p e 



in,- 



J' m {mid) 



K' (m e a) 



(3) 



where v A( > 
given by 



' (co 2 - k 2 v 2 Ai ) J m (ma) ' (co 2 - k 2 v 2 Ae ) K m ( m e a) ' 
Bo/ y/HPi, e is the internal/external Alfven speed, and the quantities m ; and m e are 

2 __ 2 (co 2 -k 2 ci)(co 2 -k 2 v 2 Ae ) 



m„ 



(4) 



1 (4 + A) {4tf - ^) ' ' (4 + vL) (co 2 - c\k 2 ) 

where c T( - e = c i( - e v A ,>/ ^J^~+v^~ is the internal/external cusp (or tube) speed. 

When an inhomogeneous transitional layer is present in the model, we cannot obtain an an- 
alytical expression for the dispersion relation unless some assumptions are made. An elegant 
method for obtaining an analytical dispersion relation is to combine the jump condi tions and the 
appro ximatio n of the thin boundar y (TB). The jump conditions were derived by 



Sakurai et al. 



(I1991a|) and iGoossens et al.1 (|1995|) for the driven problem and by Tirr v & Goossensl (I1996|) for 
the eigenvalue problem. This method was used, f or example, by bakurai et al. ( 1991b ) for deter- 
mining the absorption of sound waves in sunspots. lGoossens et al.l (119921) for determin i ng sur face 
eigenmodes in incompressible and compressible plasmas, and IVan Doorsselaere et al.l (|2004f) for 
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kink eigenmodes in pressureless coronal loops. This method is valid if the inhomogeneous length- 
scale within the reson ant layer is sufficiently small in comparison with the tube radius. Apart from 



the Alfven resonance. iSakurai et al.l(|1991af) also provided jump conditions for the slow resonance, 
which can be applied to the present situation. Hence, the TB formalism allows us to assess the 
particular contribution of the slow resonance to the kink mode damping. 

In our model the magnetic field is straight and constant so that the variations of the local 
Alfven frequency and the local cusp (or slow) frequency are only due to the variation of density. 
In that case the jump conditions at the Alfven resonance point, namely r = r A , can be written as 



m 2 /r 2 A 
ojI \d r p \j 



[£.] = -in 2 - A pr, [p r ] =0, at r = r A , (5) 



where co A and |<9 r pol A are the Alfven frequency and the modulus of the radial derivative of the 
density profile, respectively, both quantities evaluated at r = r A . The respective jump conditions at 
the slow resonance point, namely r s , are 



c 2 



oj 2 \d r po\ s \c 2 s + v 2 A 



[&] = -in ' ■ p r , [pr] = 0, at r = r„ (6) 



where to c and \d r po\ s are the cusp frequency and the modulus of the radial derivative of the density 
profile, respectively, both quantifies evaluated a, r = r„ while the fae.or (^) = (-L) . 0.034 
is a constant everywhere in the equilibrium, with y = 5/3 the adiabatic index. The jump condi- 
tions © are independent of the azimuthal wavenumber, m, but they depend on the longitudinal 
wavenumber, k z . So, for k z = 0, jump conditions © become [g r ] = [pr] = 0, meaning that there is 
no slow resonance in such a case. 

Next, we use jump conditions © and © to obtain a correction to the dispersion relation © 
due to both resonances in the TB approximation, 

: _m 2 lr\ PiPe ^k 2 I c] \ 2 p iPe 
>\ \d r p \ A w^ c 2 + v 2^ 



D m (co, k z ) = -in 2 A - in^ -f^-. (7) 

co{ \o r po\ A cof.\c 2 + v{J \d r po\ s 



The first term on the right-hand side of equation © corresponds to the effect of the Alfven reso- 
nance, while the second term is present due to the slow resonance. 



2.1.1. Expressions in the thin tube approximation 



Further analytical progress can be made by combini ng the thin boundary (T B) and the thin 
tube (TT) approximations. This was done extensively by iGoossens et al.l (|1992r) . We perform a 



-6- 



asymptotic expansion of the Bessel functions present in the dispersion relation © by considering 
the long-wavelength limit, i.e., k z a <sc 1, and only keep the lowest order, nonzero term of the 
expansion. Thus, the dispersion relation © becomes, 



m/a 



+ Pr 



m/a 



+ inpiPe 



2 I 2 

mr/ri 



1 



\drPo\/ 



+ 



i 



\d r p \ s 



= 0. (8) 



Now, we write the frequency as oj = oj r + io)\. It follows from the resonant conditions that 
CU A = CO c — 0Jr. The position of the Alfven, r A , and slow, r s , resonance points can be computed 
by equating the real part of the kink mode frequency to the local Alfven and slow frequencies, 
respectively. By this procedure, we obtain 



I 

a + - arc sin 
n 



Pi ~ Pe 



Pi +p e (A (fii+Pe) * 



(9) 



for the Alfven resonance point, and 



I 

r s = a + — arcsin 
n 



Pi ~ Pe 



2« 



Pi 



Pi+Pe 0)1 (Pi+Pe) 



(10) 



for the slow resonance point. In order to derive equation (flOl) . we have used the approximation 
u? c « c^kf, which is valid for c 2 s «; v\. Note that we need the value of cor to determine r A and r s . 
Expressions for |3 r p | A and \d r p \ s axe obtained from the density profile (eq. 0]|), 

(Pi ~Pe 



\d r po\i 



Pi-pe\X .„ . 

— - — I - cos or A , \d r po\ s 



I 



n 

-cos a s 
2 



(ID 



with a A = n(r A - a) /I and a s = n(r s - a) /I. We insert these expressions in equation ® and 
neglect terms with i.e., low-damping situation. Then we obtain an expression for the ratio 
cdr/oji after some algebraic manipulations. Since the oscillatory period, P, and the damping time, 
td, are related to the frequency as follows 



0)1 



(12) 



it is then straight-forward to give an expression for Tu/P, 



P {II a) \Pi - p e 



III 



cos a A 



+ 



(k z af 



in 



cos a. 



(13) 



where T is a numerical factor that depends on the density profile (T = 2/n in our case). As 
in previous expressions, the term with k, corresponds to the contribution of the slow resonance. 
If this term is dropped and m = 1 and cosa^ = 1, equation (TT3T ) is equivalent to equation (3) 
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of lArregui et al.1 (12008b. which only takes the Alfven resonance int o account and was first ob- 
tained by iGoossens et al.1 (| 19920 (see also equivalent expressions in iRuderman & Roberts! 12002 : 



Goossens et al 



2002) 



Next, we assume r A 
terms in equation (fT3l) is 



r s ~ a for simplicity, so cos a a ~ cos a s « 1 . The ratio of the two 



(t d )a _ (M) 
(t d ),. " 



m- 



(14) 



where (td)a and CnA are the damping times exclusively due to the Alfven and slow resonances, 
respectively. To make a simple calculation we note that the observed wavelengths of promi- 
nence oscillations correspond to 10~ 3 < k z a < 10 l . For m = 1 and k z a = 10~ 2 we obtain 
(Jd)a I ( t d). v ~ 10~ 7 , meaning that the Alfven resonance is much more efficient than the slow res- 
onance for damping the kink mode. Note that even in the extreme case of a very large 8, i.e., 

c 2 — » oo and so ( j^rj —* 1» ( t d)a / ( t d) iS ^ 1 for typical values of k z a. By means of this simple 
calculation, we can anticipate that the slow resonance will be irrelevant for the kink mode damping. 



2.2. Numerical computations 

In addition to the analytical approx imations, we also numerically solve the full eigenvalue 



problem by means of the PDE2D code (ISewelll 120051) . We consider the resistive version of the 



linear MHD equations. The magnetic Reynolds number, R m , in the corona is considered to be 
around 10 12 , but using this value requires taking an enormous number of grid points in the numer- 
ical computations. We therefore use a smaller value of R m , and consequently a smaller number of 
grid points, but we make sure that the resonant plateau is reached and the damping time becomes 
independent of the magnetic diffusivity. We use a nonuniform grid with a large density of grid 
points within the inhomogeneous layer in order to correctly describe the small spatial scales that 
develop due to both resonances. In the next Section, these numerical results are compared with the 
analytical TB approximations. 



3. Results 

We focus our investigation on the damping of the kink mode (m = 1). First, we fix the 
longitudinal wavenumber to k z a = 10~ 2 and consider l/a = 0.2. We numerically compute (see 
Fig. [2} the eigenf unctions: the velocity perturbation (v r , v v , v z ), the magnetic field perturbation (5 lr , 
B i(p , B iz ), the density perturbation (p{), and the gas pressure perturbation (pi). The behavior of v v 
and Bi r is similar to that of B\ v and v r , respectively, and hence they are not displayed in Figure [21 
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Fig. 2. — Modulus in arbitrary units of kink mode perturbations: a) v r , b) B iifi , c) p\, d) v z , e) B lz , 
and /) pi, as functions of r/a for l/a = 0.2 and k z a = 10~ 2 . The small panel in each graphic 
corresponds to an enhancement of the eigenfunction within the transitional layer. 

Because of the resonances the perturbations show large peaks in the transitional layer. The small 
panels of Figure [2] display the behavior of the eigenfunctions within the inhomogeneous layer, 
which allows us to ascertain the position of the peaks in more detail. The peaks of the perturbations 
v r , Bi v , and p x are related to the Alfven resonance, while the perturbations v,, B lz , and p x are more 
affected by the slow resonance and their peaks appear at a different position. Moreover, we see that 
peaks related to the Alfven resonance are wider than those related to the slow resonance, meaning 
that the slow resonance produces smaller spatial scales within the resonant layer. Considering the 
numerical value of co R , we get r A /a « 1 and r ja « 1.08 from equations © and (flOl ), respectively. 
Both values are in a good agreement with the position of peaks in Figure O 

Next, we plot in Figure[3]the ratio of the damping time to the period, td /P, as a function of k z a 
corresponding to the kink mode for l/a = 0.2. In this Figure, we compare the numerical results with 
those obtained from the TB approximation. At first sight, we see that the slow resonance (dashed 
line in Fig.© is much less efficient than the Alfven resonance (symbols) in damping the kink mode. 
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Fig. 3. — Ratio of the damping time to the period, td / P , as a function of the dimensionless 
wavenumber, k z a, corresponding to the kink mode for l/a = 0.2. The solid line is the complete nu- 
merical solution obtained with the PDE2D code. The symbols and the dashed line are the results of 
the TB approximation corresponding to the Alfven and slow resonances (the effect of the first and 
second terms on the right-hand side of equation Q, respectively). The shaded zone corresponds 
to the range of typically observed wavelengths in prominence oscillations. 



The individual contribution of each resonance in the TB approximation has been determined by 
solving equation CD) and only taking the term on the right-hand side related to the desired resonance 
into account. For the wavenumbers relevant to prominence oscillations, 10~ 3 < k z a < 10" 1 , the 
value of td /P related to the slow resonance is between 4 and 8 orders of magnitude larger than the 
ratio obtained by the Alfven resonance. The result of the thin tube limit (eq. lfT4lO agrees with this. 
On the other hand, the complet e numerical so l ution ( solid line) is close to the result for the Alfven 
resonance. In agreement with lArregui et al.l (2008b . we obtain Tu/P ~ 3 in the relevant range 
of k z a. As was stated by lArregui et al.l (|2008l Fig. 2d), the discrepancy between the numerical 
result and the TB approximation increases with l/a, the difference being around 20% for l/a = 1. 
However, for small, realistic values of l/a this discrepancy is less important and the TB approach 
is a very good approximation to the numerical result. For short wavelengths (k z a > 10°) the value 
of Tu/P increases and the efficiency of the Alfven resonance as a damping mechanism decreases. 
For k z a « 10 2 both the slow and the Alfven resonances produce similar (and inefficient) damping 
times. 
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4. Conclusion 



In this Letter, we have studied the kink mode damping in filament threads due to resonant 
coupling to slow continuum modes. As far as we know, this is the first time that this phenomenon 
is studied in the context of solar prominences. By considering the thin boundary approximation, 
we have found that, contrary to the Alfven resonance, the slow resonance is very inefficient in 
damping the kink mode for typical prominence-like conditions. This conclusion also holds for 
fluting modes (m > 2). The very small damping due to the slow re sonance is comparable to that 
due to thermal effects studied in previous papers (ISoler et al.ll2008|) . Therefore, we conclude that 
the effect of the slow resonance is not relevant to the damping of transverse thread oscillations, 
which is probably governed by the Alfven resonance. 
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